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We report, for the hrst time, the observation of sub- wavelength coherent image of a pure phase 
object with thermal light,which represents an accurate Fourier transform. We demonstrate that 
ghost-imaging scheme (GI) retrieves amplitude transmittance knowledge of objects rather than the 
transmitted intensities as the HBT-type imaging scheme does. 
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I. INTRODUCTION 

In many imaging circumstances, phase information 
about objects plays a role as well as or even more impor- 
tant than intensity does, for example, when the objects 
are pure-phased, that is, highly transparent and absorb 
little light, imaging can not be simply realized by the 
transmitted or reflected intensity information of thermal 
lights. Although phase distribution about an object can 
be retrieved from its Fourier-transform diffraction pat- 
tern was firstly proposed by Sayre [l| and dedicated ef- 
forts described in the works like [2] demonstrated and de- 
veloped the techniques, the efforts seems to be in vain if 
diffraction imaging applications were in hard x-ray, 7-ray, 
or other wavelengths where no effective lens or/and no 
coherent source is available. Recent works @,H| reported 
a new version of the landmark Hanbury Brown and 
Twiss (HBT) experiment @ and gave lens-less Fourier- 
transform pattern of a Youngs double slit with thermal 
sources, but the phase information is not yet mentioned 
because the object they used was amplitude-only. In fact, 
as we will discuss latterly, the classical HBT-type imag- 
ing, features that the joint detection plane in the optical 
path passes through the object, will be invalid for re- 
trieving phase knowledge about an object. Since middle 
years of last decades ^ghost imaging (GI) has been en- 
thusiastically studiedji 0, i, S 53, El, Gj, El Ell , here, 
the reason for the term ghost used is that the image of 
an object, diffractive or geometrical, would appear as a 
function of the position in the path that actually never 
pass the object, and this unique feature is regarded as 
a key difference from classical HBT-type imaging. Al- 
though whether the entangled beams was a prerequisite 
once have been hotly debated, @, H, U, Ell EH it is gener- 
ally accepted now that classical thermally emitted light 
can be used for GI and quantum entangled beams is not 
a prerequisite. In this letter, we report, for first time, 
the lensless retrieval of sub-wave length coherent imag- 
ing of a pure-phase object, which represents an accu- 
rate Fourier transform, within fresnel diffraction range 
by using GI scheme. Also, we experimentally demon- 
strate that ghost-imaging scheme (GI) retrieves complex 
amplitude transmittance of objects rather than the trans- 
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FIG. 1: Setup schemas for experiment.^ (2:, Xk) refers to the 
impulse response function of both optical paths formed by 
beam splitter from source plane X to detection plane Xt;k — 
1,2. 



mitted intensities as the HBT-type imaging scheme does. 

We first go over the theory skeleton for two-photon 
joint-detection and its related physics behind the lens- 
less Fourier transform imaging in Section [TTJ then, in Sec- 
tion [Till the experiment setup, including preparation for 
pure-phase object and pseudothermal light source and ar- 
rangement for instant intensity detecting and recording 
were briefly introduced, after that, the experiment re- 
sults were shown in Section IIVI by comparing with both 
theoretical anticipation and Fourier transform imaging 
performed by 2f system. The conclusions and discus- 
sions were finally in Section fVl 



II. THE THEORETICAL BASIS 
A. Joint detection 

FigQ] shows the setup for experiment. Intensity infor- 
mation is recorded at plane X\ and X2 in two different 
optical paths formed by a 5/5 none-polarized beam split- 
ter (BS). hi^{x, £1,2) refers to the impulse response func- 
tions from thermal light plane X to intensity detecting 
plane X± t 2, where, x and X\^. stand for positions on plane 
X and X\p.i respectively. 

The physics behind the joint detection for GI scheme 
in plane X\ and plane X2, as FigfT] shows, can be ex- 
plained in simplicity as follows: The two-photon ampli- 
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tude described by state vector \A) can be decomposed 
into weighted sum of the following three normalized ba- 
sis states: \a) = \vac, j), two-photons both reflected by 
BS; \(3) — \m, n), one photon reflected by the BS and the 
other transmitted the BS; and I7) — \i, vac), two-photons 
both transmitted the BS. i.e. \A) = ^\a) + -j=e t6l \P) + 

^e l62 \j). In the equation, 61(2) is the phase of complex 
weight for state |/?)(|7)) relative to state |a). The expres- 
sion of \a) and I7) can be expressed as: 



and 



|7> 



|wc)i|j) 5 



\th \vac) 2 , 



(1) 



(2) 



but as for the characters of identical bosons, the two- 
photon state of 1/3) must be expanded in this way: 
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-j=(\m)x\n) 2 + \n)i\m) 2 ). 



(3) 



Among the Eq.Q, the subscripts 1 and 2 of state vec- 
tor |) refer to joint detecting points in plane X\ and X 2 
respectively, and m, n stand for two undistinguishable 
photons. 

If we define E^ (ti^ 2 x\. 2 ) as the positive-frequency 
and negative-frequency components of the field at 
time-spatial point t\X\ and t 2 x 2 , suppose time t\ < 
t 2 , the joint detection probability per unit {time) 2 
that one photon is recorded at x\ at time t\ 
and another at x 2 at time t 2 , say, the square 
module of two-photon amplitude ip(tixi,t 2 x 2 ) corre- 
sponding to |^4) has been described by second-order 
Glauber correlation function [llj: G 2 (tiXi,t 2 x 2 ) = 
Tr{pE^-)(t 1 x 1 )E^-\t 2 x 2 )E^\t 2 x 2 )E^-){t 1 x 1 )}, where 
the density operator is defined as the average outer prod- 
uct of state \A): p = {\A){A\} av , i.e. , 

P=\\a){a\+\\m\ + \m~l\- (4) 

Subsisting Eq.([I])-Eq.(|4])mto second-order Glauber cor- 
relation function, and note that E^ {ti. 2 x\ i2 )\vac) = 
{vac\E^\ti^ 2 x\^ 2 ) — 0, we find that only state con- 
tributes to the joint detection as: 

G {2) {t lXl ,t 2 x 2 ) <x {(5\E^-\t x x 1 )E^{t 2 x 2 ) 

xE( + \t 2 x 2 )E(+\tiXi)\0). (5) 

On the other hand, if joint detection is only performed 
in the same plane Xi in case of HBT-type imaging, which 
we'll discuss latterly, we can regard the two-photon am- 
plitude described by state vector \A) as a pure state only 
concerning with two photons both transmitted the beam 
splitter, and thus the density operator of state \ A) comes 
to: p = \A)(A\, i.e. and joint detection probability re- 
duces to: 

G^ihx^hx^) = {A\E^\t 1 x 1 )E^-\t 2 x' 1 ) 

x#)(t 2 ^)£ (+) (*i^)l^>, (6) 



where, x\ and x[ stand for two different points in plane 
X v 



B. The relation between joint detection and 
intensity correlations 

Eq.([5]) and Eq.© states the fact that whether the 
joint detection performed in both planes X\ and X 2 , 
or in plane X\ only, the joint detection probability 
would be presented as average production of normally 
ordered positive-frequency and negative-frequency op- 
erator. Quantum detection theory shows that right 
side of Eq.Q and Eq.Q is proportional to the nor- 
mally ordered intensity correlation function of the optical 

fieldfUElM.e., 



and 



xi,t 2 x 2 ) oc (h(tiXi)I 2 (t 2 x 2 )) , (7) 



(tin, t 2 o;i) « (7i(tia;i)/2(t2a:i)). (8) 



The similarity of Eq.([7])and Eq.© means that square 
module of two-photon amplitude, can can be measured 
by correlation of intensities no matter whether joint de- 
tection is performed in plane X\ and X 2 on two different 
optical path formed by beam splitter or only in plane X\ 
on the optical path transmitted the beam splitter. 



C. Measuring the modulation to two-photon 
amplitude 

Previous worksjl^, HH have stated that if one assumes 
the time window At = t 2 — t\ = in the experiment, and 
defines A7( 2 )(n,r 2 ) = {h(n)I 2 (r 2 )) - (Ii(n)){I 2 (r 2 )), 
the impulse response functions for the two optical paths 
in FigtH ^1,2(^,^1,2) would be embedded into: 

AI^^xi, x 2 ) oc I / hl(x' ,xi)h 2 (x 7 x 2 ) 



x 



< E*(x')E(x) > dx'dx\ 2 , 



(9) 



with the assumption that the thermal state of light is 
characterized by a Gaussian field statistics [14j, among 
them x and x' denote position(s) at thermal source plane 
X in Fig[l] According to Collinss integral formula (l5j. 
there exists: 



hi(x, x\ 



dxo- 
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and 



h 2 (x,x 2 ) 



-ikd 
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(10) 



(11) 



3 




FIG. 2: The pure-phase object was prepared by etching two 
grooves with width of 150/^m and separating them by a 150/im 
un-etched area on a piece of 0.75mm x 10mm square quartz 
glass (JGS1). The other two un-etched areas with width of 
150/im are left symmetrically. The depth of two grooves is 
arranged to be 0.47/im. 

where, t(xo) is the complex amplitude transmittance of 
objects placed on plane Xq. As for GI schemed imag- 
ing, joint detection was performed on both optical paths. 
Under the condition of d = d\ + d 2 , substituting Eq. (fl"0|) 
and Eq.([TTJ) into Eq.®, we get: 

AI^( XuX2 )^\Ht[ 2 -^^}}\ 2 , (12) 

Among them T{ } denotes Fourier transform. 

In the other case, if joint detection is only performed 
at plane X\ on the optical path which passes through 
the object, then Eq. flTj)) will stand for both of two im- 
pulse response functions in Eq.(9). In this situation, the 
equation becomes: 

A/W^^Docl^lt^-^l 2 }! 2 , (13) 

\d 2 

where, X\ and x\ stand for position(s) at plane Xi . 

Based on the theory above, the following experiment 
was carried out in the regime of large number photons to 
illustrate the behaviors of two-photon interference. Now 
we can easily envision the physics behind the following 
experiment: The modulation to the two-photon ampli- 
tude, %l)(x\,x 2 ) and ip(xi,x[), by the objects was mea- 
sured by correlation function of intensity fluctuations. 



III. THE EXPERIMENTAL SETUP 

The experiment is set out with the use of pure phase 
object , which was prepared by etching two grooves with 
width of 150/Uto and separating them by a 150/xm un- 
etched area on a piece of 0.75mm x 10mm square quartz 
glass (JGS1), as Figj2] shows. The other two un-etched 
areas with width of 150/xm are left symmetrically. Since 
we are using Nd: YAG laser with the wavelength of 
A = 0.532//m , the depth of two grooves is arranged to be 
A/2(n - 1) = 0.532/2(1.57 - 1) = 0.47^m to form phase 
differences of A0 = n from those un-etched area. The 
thermal light source at plane X is simulated by project- 
ing a beam of laser pulse (frequency doubled Nd: YAG 



impulse laser, A = 0.532/im) onto a slow-round ground 
glass disk, the diffracted dynamic speckle field features 
the full statistical character of real thermal light but with 
much longer coherent time to allow detecting its inten- 
sity fluctuations. Away from thermal light source plane 
X 40mm, we place a beam splitter to form two optical 
paths for imaging system referring to Fig[T] shows. The 
pure phase object t(xo) is placed on plane Xq away from 
source plane X by d\ = 60mm. Charge Coupled Device 
camera CCD-1 is placed at plane X\ in the optical path 
after the object by d 2 = 75mm. Correspondingly, we 
place CCD-2 at plane X 2 on the other optical path, ar- 
range the distance from the thermal light source X to X 2 
to be d = 135mm. The parameters are thus arranged to 
meet with the condition d = d\ + d 2l which was required 
by the theoretical prediction in Ref . [i~3| . The light spot 
diameter on ground glass is chosen as a = 3mm in or- 
der to ensure linear dimensions of the coherence area(20l| 
of the pseudo-thermal light field across the object plane 
D w Xdx/a = 0.532x60/3= 10.64^m to be much smaller 
than the pure object's feature size 150/xm. The exposure 
time of two CCD cameras is set to 1ms. Totally about 
10, 000 frames of independent two-dimensional instant in- 
tensity distribution data of speckle fields in the planes of 
X\ and X 2 , say, ensembles of Ii(xi) and I 2 (x 2 ), are 
recorded to prepare for correlation operation. The laser 
pulse shooting, data acquisition and their recording pro- 
cess are synchronized and accomplished by computers. 

IV. EXPERIMENTAL RESULTS 

Choosing symmetric positions x 2 and x\ — ~x 2 to 
calculate A/( 2 )(xi,X2) with preserved data I\(x\) and 
I 2 (x 2 ), we find what we obtain (Fig[3ji) shares the same 
pattern in Fig. ([3b), but have double coordinate scales 
just as equation: 

AI^x 2 ,-x 2 )^Ht( 2n[X \-^ X2 h I 2 

= \Ht(^)}\\ (W) 

which derived from in Eq. (|12[) . anticipated. The re- 
sult (Fig[3£i) shows a sub- wavelength interference pattern 
equals to the Fraunhoufer diffraction pattern(Figj3jD) of 
the same pure-phase object accomplished by a 2/ sys- 
tem(with A = 0.532^m and / = 75mm), but with half of 
the wavelength of coherent illumination. In Fig. red 
curve represents profile of two dimensional experimen- 
tal result, and blue curve presents the standard Fourier 
transform of the pure-phase object's complex amplitude 
transmittance under the sub-wavelength condition ac- 
cording to Eq. (fl"4]) . Its analytical form is given by 

AI 2 (x 2 , ~x 2 ) = sinc 2 ( a x 2 )[l - 2 

XsubO-2 

2ira 4ira 9 
x cos(- —x 2 ) +4cos(- -7-2:2)] , (15) 

A S m6«2 ^subU 2 
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FIG. 3: a. Retrieved two dimensional patterns of AI(x2, —X2) 
and its profile (red curve) by using GI schemed setup; b. The 
two dimensional Fraunhoufer diffraction pattern and its pro- 
file of the pure phase object obtained by coherent 2/ system 
with / = 75mm and A = 0.532/xm. 



in which a = 150/xm, dn — 75mm, and X su b = | 32 /im. 
From FiglHK, we can see the experiment result fit the 
theoretical Fourier transform Eq. (|15|) quite well. 

According to Ref.[5(, the optical path that passes the 
object in our experimental setup is compatible with the 
classical HBT-type imaging system if joint detection is 
accomplished only on the arm passes the object. So, 
the 10,000 frames of previously preserved intensity data 
recorded on plane X\, i.e., ensembles of I{xi) can also be 
utilized for HBT type of imaging. In this situation, the 
two points x\ and x' x for joint detecting are chosen only at 
plane Xi symmetrically with respect to the center of it, 
i.e., x\ = —x[. The function AI 2 (xi,— x\) contains no 
phase information about the object in this time, as two- 
dimensional pattern and its profile (red curve) in Fig0] 
shows. Eq. (fl~3|) gives a straightforward explanation to the 
experimental result: The diffraction pattern we obtained 
was the Fourier transform of the square module of the 
transmitted function of the object. The blue curve in 
FigH] describes the Fourier transform of |£(xo)| 2 , with 
the analytic form of 



AI 2 {x\, —xi) = sinc 2 ( 
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-xi) 



(16) 



which also fits the experimental result very well. 

Thus we find the essential differences of our GI 
schemed sub- wave length coherent imaging from recently 
reports about sub-wavelength interference with thermal 
light 0, H[ is what GI type of imaging system in our 
experiment setup retrieved was the complex amplitude 
transmittance knowledge of the object rather than the 
transmitted intensity as the HBT schemed imaging does. 



FIG. 4: Retrieved two dimensional pattern of AI(xi,—xi) 
and its profile by using HBT-type of setup, recovers a diffrac- 
tion pattern of the limited aperture of the object and retrieves 
no information of the phase object. 



CONCLUSION AND DISCUSSION 



The lensless scheme proposes the potential application 
in hard x-ray, 7-ray, or other wavelengths where no ef- 
fective lens or /and no coherent source is available, and 
the similar idea have also been reported brilliantly, but 
in spatial domain[2l|. On the other hand, obtains of sub- 
wavelength interference pattern suggest that diffraction 
limit can be broken through by using two-photon absorp- 
tion (TPA) media if we find ways to fold the symmetric 
planes of x\ and xn into the same one22[. This fea- 
ture appeals to quantum lithography which was widely 
discussed 231 . 



Apart from the unique sub-wavelength feature of the 
experiment result, attention shall also be paid again to 
the similarity and the difference between GI schemed and 
the classical HBT-type of imaging: They both retrieve 
diffraction patterns by correlation function whereas GI 
recovers knowledge of complex amplitude transmittance 
about objects rather than the transmitted intensities as 
HBT-type imaging does. Similar conclussion have also 
been reported 24( , but with use of lens in both optical 
paths. In fact our scheme suggests a new way to per- 
form lensless Fourier transform within fresnel diffraction 
region. 

Unlike the other ghost imaging and ghost diffraction 
experiments, a pulsed thermal-like source is used instead 
of a continuous one, the intensity correlation can be mea- 
sured even the exposure time of CCD camera is much 
longer than pulse width of the source. This unique im- 
pulse feature of our thermal light source enlightens us 
on the issue for recording intensity fluctuation by slow 
detector. This way suggests achieving to record inten- 
sity fluctuating much faster than the respond speed of 
the detecting system. As for hard a;-ray imaging, there is 
a potential applicability to record the femto-second fluc- 
tuations of intensity by using of detecting system with 
response speed in nanoseconds. 
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